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The total momentum of N interacting bosons or fermions in a cube equipped with periodic boundary 
conditions is a conserved quantity. Its eigenvalues follow a probability distribution, determined by 
the thermal equilibrium state. While in non-interacting systems the distribution is normal with 
variance ~ N, interaction couples the single-particle momenta, so that the distribution of their 
sum is unpredictable, except for some implications of Galilean invariance. First, we present these 
implications which are strong in ID, moderately strong in 2D, and weak in 3D. Then, we speculate 
about the possible form of the distribution in fluids, crystals, and superfluids. The existence of 
phonons suggests that the total momentum can remain finite when N —>■ oo. We argue that in 
fluids the finite momenta distribute continuously, but their integrated probability is smaller than 1, 
because the momentum can also tend to infinity with N. In the fluid-crystal transition we expect 
that the total momentum becomes finite with full probability and distributed over a lattice, and that 
in the fluid-superfluid transition a delta peak appears only at zero total momentum. Based on this 
picture, we discuss the superfluid flow in both the frictionless and the dissipative cases, and derive 
a temperature-dependent critical velocity. Finally, we show that Landau’s criterion for excitations 
in moving superfluids is an in some cases correct result of an erroneous derivation. 
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Our concern in this paper is to understand the role of 
the total momentum in the description of fluids, crystals 
and superfluids. The idea that such a role can exist at 
all is based on the experimental fact that in interacting 
quantum systems there are excited states that involve 
a large number of particles and carry a finite total mo¬ 
mentum. Such collective modes exist at all energies, in 
all the thermodynamic phases; we think of compression 
waves or lattice vibrations [1]. This fact is in a striking 
contrast with the naive expectation that, because the to¬ 
tal momentum is the sum of N single-particle momenta, 
it should be normally distributed with a variance ^ N. 
Deviation from this can be a joint effect of interaction 
and quantum mechanics. In interacting quantum sys¬ 
tems the single-particle momenta are not independent: 
because they do not commute with the interaction, free 
energy minimization couples them. Therefore, the distri¬ 
bution of their sum is largely unpredictable, although we 
shall see that Galilean invariance does have implications 
for it. In classical systems the single-particle momenta 
are identically distributed independent random variables, 
so their sum is normally distributed with variance ^ Af, 
and similar holds in non-interacting quantum systems. 

Once the total momentum can be either finite or tend¬ 
ing to infinity with N, one can speculate about the 
changes that can take place in phase transitions. Two 
characteristic changes are to be expected when the tem¬ 
perature T decreases. First, the continuous distribution 
for the finite momentum, certainly existing at high T, 
can partly or entirely be replaced by a discrete one. Sec¬ 
ond, the total momentum of order '/N, also present at 
high T, can disappear. We shall argue that both occur 
in the fluid-solid transition, while in the fluid-superfluid 
transition only a delta peak appears at zero total mo¬ 
mentum. We do not expect a qualitative change in the 


vapor-liquid transition. 

A crucial point that must be stressed is that the to¬ 
tal momentum is a property of the entire system. Saying 
that in He H below the A-point it can be zero with proba¬ 
bility p and nonzero with probability 1—p will mean that 
the equilibrium state is a convex combination of a super¬ 
fluid and a normal fluid state with respective weights p 
and 1—p. This is the least intuitive aspect of our pro¬ 
posal: that a physical meaning is attributed not to parts 
of the system (e.g. to Nq < N particles forming the 
Bose condensate) but to parts of the density matrix. A 
detailed account of this work is given elsewhere [2]. 

Total momentum and Galilean invariance.-As in our ear¬ 
lier paper [4], we consider the energy operator 

I ^ 

of N interacting particles in a d-dimensional cube A of 
side length L, defined with periodic boundary conditions, 
which admit Galilean transformation in a bounded do¬ 
main. C7 a is the periodized potential energy containing 
translation invariant pair and possibly many-body in¬ 
teractions, which can include hard-core interactions but 
no external field. In such a case the total momentum 
P = Pj = —ifiJ2j is conserved and the eigen¬ 

values and eigenvectors of H can be labeled by the eigen¬ 
values hQ of P: PV’Q.n = flQ^Q.n, 
where Eq ^ < EQ^n+i [n > 0) and 

Q e A* = {(2ni7r/L)ei H-h {2ndTT/L)ed} 

( 2 ) 

'i/’o.o and Eq q are the ground state and its energy. A 
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iV^-point subset of A*, that we call irreducible, 

^bred = {q G A* : -ttN/L < qi< ttN/L, i = 1,. .. ,d}, 

( 3 ) 

plays a distinguished role: any Q G A* can uniquely be 
written as q + A^k, where q G At and k G A*, and 
all the eigenvalues and eigenvectors are obtained by a 
Galilean boost from those belonging to [4], 

Eq+Nk.n = Eq^n + {h^/2m)[N'k^ + 2k • q], 
V'q+wk.Ti = q,kGA*,n>0, (4) 

where x = N~^ '^J=i center of mass. With j3 = 

and Z = X]q the density matrix 

is ^"^EQE„e“^^‘^’”IV’Q.n)(V'Q,nl- The sum 






-PEct,. 


( 5 ) 


n—0 


is interpreted as the probability that P = ?iQ in thermal 
equilibrium. Taking the thermodynamic limit at a fixed 
density p = N/L'^, one can prove the following [2]. 

(i) In ID, with \p = yj2TTl3W'/m, 

lim Prob I < xl = [ 

JV^oo 1^ j TT Jo 

(ii) In 2D, for any J > 0 integer, 


e dy. (6) 


lim Prob 

N —¥GO 




^ —const. X 


(iii) For d > 3, 


( 7 ) 


lim ProbilQd < 3nN/L} = I 

AT—>00 


( 8 ) 


[but we believe that |Q|/-\/]V —>• 00 has a vanishing 
asymptotic probability, just as for d = I, 2]. 

One can add that for non-interacting systems a nor¬ 
mal distribution of variance ~ — Nq [Aq is the num¬ 

ber of particles with zero momentum] is expected in any 
dimension. The specificity of ID is that the irreducible 
wave numbers are in the interval (—Trp, 7rp] independently 
of the system size. Because the interaction fully re¬ 
veals itself in the irreducible ensemble, this gives a hint 
to why, on the scale Q ~ Vn, do any ID system re¬ 
semble a non-interacting one [3]. On the other hand, 
maxqgA*^^^ IQij oc \/W in 2D and oc A^/^ in 3D; for 
this reason one cannot expect that the analog of Eq. (6) 
would hold for interacting systems in 2D and 3D. 

Density matrix reduced to the center of mass.-P is canon¬ 
ically conjugate to x. A consequence is that the depen¬ 
dence on X can be separated in V'Q.n: with the notations 
X = (xi,.. .,xjv), X' = (x^,... ,x'^), x' = Xj - xi. 


V^Q.„(X) = 


_ piQ-x^-(*/W)Q-EL 


-^"2^q,„(0,X'). (9) 


We can reduce the density matrix to the center of mass 
through a partial trace over X' [5], with the result 

(xjpc.m.ly) = P (QIPc.m.lQ') = 

qga* 

( 10 ) 

so Ni/q are the eigenvalues of Pc.m.) with corresponding 
eigenfunctions (x|Q) = L~ie^^'^. 

Equations (4)-(I0) are valid at all T and p, for bosons 
and fermions, with or without interaction. From now on, 
we deal only with bosons. We recall that the one-particle 
reduced density matrix is 

(x|pi|y)=p (Q|pi|Q') = 

QGA* 

( 11 ) 

where nq = (Aq)/A and Aq is the occupation num¬ 
ber operator for the single-particle state |Q). Because 
XqeA* = 1 just as XqeA* = 1, "Uq can be 
interpreted as the probability that a particle is in the 
one-particle state |Q). Thus, there is a perfect struc¬ 
tural analogy between pc.m. and pi. They have the same 
eigenvectors |Q), the corresponding eigenvalues are prob¬ 
abilities multiplied by A, the expressions (10) and (II) 
are nonnegative, therefore vq and nq are positive defi¬ 
nite functions on A*: the n x n matrices [nq._q^]ij and 
[nQ^-Qj]ij are positive semidefinite for any n and any 
Qi,..., Q„ in A* [2]. In particular, 

^0 > vq, no > nq any Q G A*. (12) 

We note that vq > nq implies that the free energy density 
can be computed in the Q = 0 ensemble [2]. 

Distribution of the total momentum in infinite space- 
Structural similarity between pc.m. and pi is contrasted 
with a fundamental difference, owing to the fact that the 
dependence of the eigenstates V’Q.n on x is separable, but 
the dependence on Xj is not. The difference is accentu¬ 
ated in the thermodynamic limit in which both nq and 
nq tend in distribution sense to a positive and positive 
definite measure on the dual space . We denote 
them by n(k) and h(k), respectively. If there is Bose- 
Einstein condensation (BEG) then h(k) = no(5(k)-|-hc(k) 
where no > 0 and he is continuous. Now f h(k) dk = I 
always holds true, otherwise the kinetic energy density 
would diverge in the thermodynamic limit [6]. So h is a 
probability measure, meaning that in infinite space, as in 
finite volumes, the single-particle momentum cannot be 
infinite. This is not true for h: in finite volumes the total 
momentum is finite with probability 1, but from Eqs. (6) 
and (7) one can infer that p<oo = / f'(k) dk = 0 in ID 
and p<oo < I in 2D for all T > 0; also, p<oo < 1 should 
hold for high T or small p in 3D. p<oo and Poo = 1 —p<oo 
are the respective probabilities that the total momentum 
is finite or infinite in infinite volume. 

p<oo = 0 means that states carrying a finite total mo¬ 
mentum, compression waves included, gradually lose all 
their statistical weight with the increasing system size. 
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This occurs at T > 0 in non-interacting systems in any 
dimension, and also in interacting systems in ID. This 
does not mean that density waves as finite-energy excita¬ 
tions do not exist. However, in any dimension the states 
with an energy gap of order 1 above the ground state can 
play a role only at T = 0: at T > 0, i?Q,n — i?o,o cx; N 
for the relevant eigenstates. 

In contradistinction to collective excitations, few- 
particle excitations involve small groups of particles 
which are nearly separated from each other. Rotons may 
be of this kind [7]. In 3D at high energies, in 2D at all en¬ 
ergies they are expected to occur with a positive density 
in a non-vanishing fraction of eigenstates, and their finite 
momenta to add up randomly to a total momentum of 
order y/N. 

Total momentum and thermodynamic phases: 
Speculations.-Now we can make our proposals about the 
connection between i> and the thermodynamic phases. 
The above two classes of excitations, the collective 
ones with a finite nonzero total momentum and the 
single- and few-particle excitations with an accumulated 
infinite total momentum in infinite space are sufficient 
to characterize the fluid phases (gases and liquids) at 
T > 0 and d > 1. Namely, D is continuous, and the 
equilibrium state S in infinite volume is of the form 
S = P<oo '5<oo +Poo Soo- Here 5<oo and Soo themselves 
are composed of states with finite or infinite momentum, 
respectively, and can be obtained via suitable limits from 
parts of the density matrix [2]. Since in 2D the second 
term is nonzero at all T > 0, a 2D system cannot be in a 
crystalline state (see below). This corresponds to known 
results about the preservation of shift-invariance [8-11]. 

Concerning crystals, physical intuition suggests that 
the random single- and few-particle motions freeze out in 
the thermodynamic limit, implying p<oo = 1- Moreover, 
the limit of (xjpc.m. |y) must be periodic in y— x according 
to some lattice L, therefore v{\T) = X^ksl* ^K'5(k — K), 
where L* is the reciprocal lattice. Because of positive 
definiteness, vq > f'K- Coherent quantum crystals [12, 
13] are curious objects interpolating between crystals and 
fluids. For a discussion see Refs. [2], [6]. 

The characterization of superfluids raises two ques¬ 
tions: what is the role of BEC, and what are the two 
fluids in Tisza’s two-fluid picture [14]. There is today a 
consensus based on analytical, numerical and experimen¬ 
tal results [15-21] that in bulk liquid helium BEC and 
the superfluid transition occur simultaneously. Yet in 2D 
there is superfluidity [22, 23] but no BEC at T > 0 [24- 
26], and bulk helium at T = 0 is 100% superfluid con¬ 
taining less than 10% of condensate. These controversial 
facts can be reconciled by assuming that both in 2D and 
in 3D the primary event at the superfluid transition tem¬ 
perature Tg is ’condensation’ of the entire system into a 
state of zero total momentum, and in 3D this triggers 
settling of individual particles at zero momentum. Thus, 
when passing Tg, the continuous distribution v picks up 
a Dirac-delta, so that it becomes f'(k) = f'o<5(k) -I- £'c(k), 
where vq = lim^^o bniL^oo Z[|q|<K > 0 and Vc is con¬ 


tinuous. An immediate implication is that the super and 
normal fluids are not parts (as the condensate) but al¬ 
ternative states of the whole system, that compose the 
equilibrium state by convex combination, 

5 = vo{T) Sg{T) + [1 - vq{T)] Sn{T). (13) 

The super and normal fluid fractions are vq and 1 — f^o, 
respectively. Sg is a state of zero total momentum, and 
Sn can further be decomposed according to whether the 
total momentum is nonzero finite, or infinite. Eor non¬ 
interacting bosons £'(k) = 0 at T > 0, so in spite of BEC, 
their state is not a superfluid. 

Frictionless /lotc.-Galilean boost generates frictionless 
flow. If the flow has a velocity v with respect to the 
container A, and mv/fi € A*, then in the comoving ref¬ 
erence frame we have [2, 4] 

Tt V^Q,n — AQ_7Ymv/?l,nV^Q,n; 

= (SQ - A^mv)'i/'Q,„. (14) 

H'' and P'^ are obtained from H and P, respectively, 
by replacing pj with — mv [27]. H'' has the same 
eigenvalues and eigenfunctions as H: the velocity boost 
only cyclicly permutes the eigenstates among the ener¬ 
gies, by adding Nmv to their momentum. The density 
matrix reduced to the center of mass will be (0]/?^ ^^^ jx) = 
e ‘fi ""'"■’’^(Ojpc.m.jx), and its square root can play the role 
of a macroscopic wave function. Eor large L and large jxj, 
(OjPc.m.jx) « pf'o- The super and normal fluid states in 
infinite volume can be constructed, and Eq. (13) becomes 

S^{T) = MT) S^iT) + [1 - i>o{T)] S:{T), (15) 

with unchanged super and normal fluid fractions. This 
is an equilibrium state whose free energy density is the 
same as that of the fluid at rest. 

Dissipative flow: Speculations.-A real flow is frictionless 
in the superfluid state but is dissipative in the normal 
fluid state. If in the superfluid state the velocity is v, 
then in the normal fluid state it will be smaller and time 
dependent, a{t)-v, a{t) < 1. The difference is due to the 
produced heat which partly is dissipated into the envi¬ 
ronment, partly raises the temperature to Tt > T in both 
states [28]. Now S^{T) is replaced by 

SLne^iTt) = MTt) SJiTt) + [1 - f>o(rt)]5“(‘)''(T,). 

( 16 ) 

This is a non-equilibrium state because the two terms 
are generated by different Hamiltonians, and 
respectively [29]. Tt can be computed from 

^ksTt - ^ksT = - a{t)% (17) 

where 77 < 1 is the efficiency, i.e., the fraction of the heat 
that raises the temperature. The dissipation continues 
until a(t) = 0. It raises Tt and drives cx{t) to zero. In the 
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limit t —>■ c» the system goes over into a non-equilibrium 
steady state at a temperature Ty, 

SLne^iTy) = MTv) 5J(r„) + [1 - Po{Ty)] 5„(T„). (18) 

From Eq. (17), Ty = T + rimv‘^/U b- As v increases, T„ — >• 
Ts, t'o(T'-u) —>■ vo(Ts) = 0, and the system ends up in the 
thermal equilibrium state iS„(Ts), at rest with respect to 
A. The critical velocity can be inferred from the equation 

= Ts which yields 

v„{T) = ^kB{Ts-T)/i7jm). (19) 

For He II, with Ts = 2.17 K, ?;cr(0) = (67/0?) m/s > 67 
m/s. 

One can qualitatively understand why friction affects 
the moving normal fluid but not the superfluid. Fric¬ 
tion (viscosity) can be described as a spatially random 
uncorrelated perturbation. The eigenstates forming the 
superfluid state must be strongly entangled and, thus, 
resistent to such type of perturbations. Friction can only 
influence small separated groups of particles and excite 
them, at the expense of the kinetic energy of the flow, 
to a higher energy. Because its effect extends to the en¬ 
tire volume, the overall increase in energy is of order N. 
When these groups relax, they emit incoherent radiation 
a part of which leaves the system, another part is reab¬ 
sorbed and heats it. The suitable separated small groups 
of particles can be found in a macroscopic number only 
in the eigenstates ipQ+Nmv/h,n = with 

IQI ^ y/N, which contribute to S^. 

Comment on Landau’s critical velocity.-We show that 
Landau’s critical velocity [30-32] was obtained through 
an erroneous argument. In Section 4 of Ref. [30], Lan¬ 
dau investigated the stability of the superfluid ground 
state of liquid helium in a capillary against low-energy 
excitations that a flow of velocity v can create by losing 
kinetic energy. His argument is based on the properties 
of €q = Eq^o — Eq^o, the energy gap to the lowest-lying 
eigenstate of momentum fiQ. In He H, eg is measured 
up to about 4A“^. Its qualitative features do not de¬ 
pend on T(< Ts): the curve starts linearly, passes over 
a maximum and exhibits the roton minimum [33-35]. If 
eg could be measured at T = 0, it would certainly show 
the same features as in the measurement [33], done at 1.1 
K. In our notations, Landau wrote down the equation 

eq+Nm^,/h = Cq -b ftq • V -k (20) 

first with Cq « cli|q| near q = 0 [Eq. (4,1)] and sec¬ 
ond, with Cq « A -b (2/r)“^?i^(q — q^.)^ [Eq. (4,3)] near 
the roton minimum q^ [36]. He considered Eq. (20) as 
describing the energy balance of the moving fluid, so, ac¬ 
cording to him, to excite the system at the expense of a 
part of the kinetic energy, 

£„ -b ?iq • V < 0 or v > for some q (21) 


must hold. In particular, rotons are excitable if v > 
fqr/ (^lli-j), and phonons are excitable iiv > c. However, 
the velocity threshold (21) has no special meaning in a 
flow experiment. If we add to both sides of Eq. (20), 
we obtain 

Aq+ATmv/fi.o = Aq,o + ftq ' V -b (22) 

the first of Eqs. (4) with n = 0 and k = mv/h G A*. 
By Eq. (22), the set {Eq_o -b ?iq • v -b ^Amv^jqgA* is 
the same as {Eq^olqeA*, with a permutation of the el¬ 
ements, which changes as v varies. Since the energy 
spectrum is unchanged, both the ground state energy 
and the free energy are independent of v. As seen from 
Eq. (14), for any v the minimum energy Eo,o is assigned 
to the lowest-lying eigen¬ 

function that describes the system moving with velocity 
V, while the energy of ^o.o is increased by iVmv^/2. They 
only exchanged energies, Landau’s flow is nondissipative. 
In flow experiments eg plays no role: when the normal 
fluid is present, the real flow is dissipative and, at any 
r: > 0, its loss of kinetic energy is of order N , highly 
above eg from |Q| = 0 to IQI ~ '/N . Thus, if eg should 
be compared with the available excitation energy to get 
the critical velocity, then this would be zero. Neverthe¬ 
less, condition (21) is relevant when bulk He H is locally 
perturbed by moving in it a tiny object (e.g. an ion) with 
velocity v. Then, the energy transfer is indeed of order 1, 
and rotons and phonons are created above the respective 
velocity thresholds. 

Final note.-At 0 < T < T^, eg is qualitatively the same 
as at T = 0, but the energy of the state at Q = 0 is in the 
order of N higher than the energy of the ground state. 
This must be the superfluid state at the given temper¬ 
ature, providing an experimental support to our claim 
that the superfluid is a state of zero total momentum. 

Summary.-We studied the distribution of the total mo¬ 
mentum of interacting quantum systems in continuous 
space. We derived a few mathematical results and com¬ 
bined them with an experimental fact, the existence of 
finite-momentum excitations, to make a bold intuitive 
jump into suggesting that the thermodynamic phases 
could be identified on the basis of the distribution of the 
total momentum in the limit of infinite volume. Different 
total momenta correspond to different states whose sum, 
weighted with their probabilities, gives the density ma¬ 
trix. Suitably devised measurements can project the sys¬ 
tem into states present in the density matrix. We made 
specific suggestions about the momentum distribution of 
fluids, crystals and superfluids, discussed the frictionless 
and the dissipative flow of the latter, derived a formula 
for the critical velocity, and pointed to a failure in the 
derivation of Landau’s criterion. 
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